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Abstract 

Based on a simple model for spin-polarized scanning tunneling spectroscopy (SP-STS) we study 
how tip magnetization and electronic structure affects the differential conductance (dl/dV) tunnel- 
ing spectrum of an Fe(OOl) surface. We take into account energy dependence of the vacuum decay 
of electron states, and tip electronic structure either using an ideal model or based on ab initio 
electronic structure calculation. In the STS approach, topographic and magnetic contributions to 
dl/dV can clearly be distinguished and analyzed separately. Our results suggest that the sensi- 
tivity of STS on a magnetic sample can be tuned and even enhanced by choosing the appropriate 
magnetic tip and bias set point, and the effect is governed by the effective spin-polarization. 
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I. INTRODUCTION 



Development of simulation tools for scanning tunneling microscopy (STM) and spec- 
troscopy (STS) are in the focus of theorists [J since the invention of STM 30 years ago. 
It has been established that the tip electronic structure plays a crucial role in measured 
differential conductance (dl/dV) tunneling spectra, see e.g. Refs. js, 4|. While a theoretical 
method has been proposed to separate tip and sample contributions to STS jsl, some recent 
research activities have focused on extracting surface local electronic properties from exper- 
imental STS data sl IgUqI, which is the convolution of tip and sample electronic structures. 
The situation is expected to be even more complicated in magnetic systems due to effective 
spin-polarization. 

STM can be made sensitive to magnetism, and spin-polarized scanning tunneling mi- 
croscopy (SP-STM) and spectroscopy (SP-STS) are nowadays the key tools for nanomag- 
3Nl2|. i.e. for studying and manipulating magnetic properties of surfaces and 



netic research 



10- 



deposited magnetic nanoclusters with atomic scale resolution |l3l-ll5|. SP-STS has recently 



been used to find inversion of spin-polarization above magnetic adatoms 



16-18 



which effect 



has been explained theoretically 19|. This spectroscopic approach turned out to be useful 



for investigating many-body effects on substrate-supported adatoms 20|, l2l|] as well. Tip ef 



fects on SP-STS 



22| and on achieving giant magnetic contrast 



23 



have also been reported. 



Recent experiments of Schouteden et al. {2^ and Heinrich et al. [25| show evidence that STS 
peaks are different above magnetic islands of opposite magnetization. However, no detailed 
explanation has been given for these observations. In our present study we shed light on 
this effect in a different setup by studying tunneling spectra of the magnetic Fe(OOl) surface. 
Our results show the importance of the effective spin-polarization on the STS spectra, which 
can be tuned by changing tip magnetization direction or bias set point for the tip. 

The paper is organized as follows: Theoretical model of SP-STS is presented in section ITTl 
Based on this, simulating differential tunneling spectra of the Fe(OOl) surface is presented 
and discussed in section IIIII We focus particularly on tip effects, and consider an ideal 
electronically fiat and maximally spin-polarized tip as well as a more realistic ferromagnetic 
Ni tip. Our results suggest that the sensitivity of tunneling spectroscopy measurements on 
magnetic surfaces can be enhanced by using the proper magnetic tip and bias set point. 
Summary of our findings is found in section IIVI 



2 



II. THEORETICAL MODEL OF SP-STS 



Based on the Transfer Hamiltonian approach, the work of Passoni and Bottani 26| re- 
ports an advanced way to incorporate results of electronic structure calculations into STS 
simulations. Our model for SP-STS is based on the spin-polarized version of the Tersoff- 



Hamann model 



27 



28 



introduced by Wortmann et al. 29|, which is adapted within the 



28 



30 



-|32[. App roximations and limitations 



28 



3l|. By simulating differential 



framework of the atom-superposition approach 
of this approach have been discussed in the literature 
conductance, we directly calculate the differential spectrum and not the numerical derivative 
of an integral spectrum, see Ref. [sl- This way we assume that dl/dV is proportional to the 
electron local density of states (LDOS), which is a reasonable approximation at low bias 
voltages. Calculating dl/dV from the tunneling current by numerical differentiation is a 



commonly used approach 



3 fly, 



26(1 . The motivation comes from experiments, where the 



surface LDOS is not known, and it is the goal to extract this quantity from measurements. 
In our SP-STS approach based on first principles electronic structure calculations we pro- 
ceed in the opposite way: We calculate surface electronic structure, define LDOS at tip 
apex position, also taking into account tip electronic structure, and then we define dl/dV 
proportional to LDOS. This way, the so-called background term does not explicitly occur 
in our dl IdV expression, which is an important correction to dl/dV at higher bias voltages. 



see Ref. 



3|, Eq.(3), second term. Assuming an electronical 



y featureless tip, this term is 
81], and, thus, it is expected to 



approximated to be proportional to the tunneling current fl, 
increase with increasing absolute value of the bias voltage. Due to the energetic position of 
the Fe(OOl) surface state peak close to the Fermi level, in the present work we neglect the 
background term and we will deal with its relevance in the future. We expect, however, that 
the inclusion of the background term does not affect our conclusions for tuning the dl/dV 
peaks by using a magnetic tip as long as staying in the low bias regime. Resolving dl/dV 
features at higher bias voltages turned out to be difficult even for using nonmagnetic tips 
and, therefore, different normalization schemes have been introduced to obtain information 



about the sample local electronic structure [3|, [6 



i9i. 



Incorporating bias dependence into our model one has to take into account energy de- 
pendence of the vacuum decay of electron states. Moreover, we show how to explicitly 
incorporate energy dependence of tip electronic structure based on the result from ab ini- 



3 



tio calculation into our bias dependent atom-superposition-based model. This means that 
different tip models and their effect on tunneling properties can be investigated. The only 
requirement for our present formalism is that we assume that electrons tunnel through one 
tip apex atom. 

In order to simulate single point differential tunneling spectra above the surface atom 
with lateral coordinates {xo,yo), first, we calculate the LDOS at the tip apex position 
RTip{xo,yo, z). The spin-mixed LDOS at a given energy E can be decomposed into a 
non-spin-polarized (TOPO), and a spin-polarized (MAGN) part as 

LDOS{xo, yo, z, E) = LDOStopo{xo, yo, z, E) + LDOSmagn{xq, yo, z, E). (1) 



Using the atom-superposition method (3l|, |32|, the two terms can be written as 

LDOSTOPo{x,,y,,z,E) = e-2''(^)l^"^(^«'^°'^)-^<^lnT(^)ng(^) (2) 

a 

LDOSMAGN{xo,yo,z,E) = A^^e-2'^(^)l^"^(^«'^«'^)"^<^lmr(i^;)mg(E)cos</.„, (3) 

a 

where the sum over a has to be carried out over all the surface atoms with position vectors 
R^, each characterized by a local spin quantization axis determined by their atomic spin 
moment direction. AE is the energy resolution for our simulated tunneling spectra, and it 
ensures that the LDOS is correctly measured in units of 1/eV. A AE value of 10~^ eV 
has been used in our calculations. The exponential factor is the transmission coefficient for 



electrons tunneling between states of atom a on the surface and the tip apex, w^ 
the vacuum decay, k, is treated within the independent-orbital approximation [27 



here k is 



28 



which means that the same decay is used for all type of orbitals, but its energy dependence 



is explicitly considered in the same fashion as in Ref. 33|. Extension of our model in 



3^ is 



the direction to incorporate orbital dependent vacuum decay following Chen's work 
planned in the future. In the present paper we propose two different ways of calculating k, 
one is inspired by the Tersoff-Hamann model, taking only surface properties into account. 



k{E) = ^^2m{<Ps + E'p-E), (4) 

where the electron's mass is m and charge — e, while (ps and Ep are the average electron 
workfunction and the Fermi energy of the sample surface, respectively. We use this energy 
dependent vacuum decay for an ideal, electronically featureless and maximally spin-polarized 
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tip model. The second expression for k is based on the Wentzel-Kramers-Brillouin (WKB) 
approximation assuming a rectangular tunnel barrier, 




2™(^^i^.^|-s). (5) 

with ipT being the local electron workfunction of the tip apex, and V is the applied bias 
voltage. This vacuum decay formula is considered for our magnetic Ni tip model. The 
quantity {(ps + ft^T + ^V) /2 + Ep — E is the energy- and bias dependent apparent barrier height 
for tunneling electrons, (paiE, V). Note that our vacuum decay formulae are asymmetric with 



respect to positive and negative bias regime jo]. The average workfunction of the sample 
surface is calculated from the local electrostatic potential on a three-dimensional fine grid, 

V{x,y,z), as 

' ~ " (6) 



max 



y x,y J 



with Nx and Ny the corresponding number of grid points, and the local workfunction of tip 
apex is obtained as 

(j)T = max {V{xo, yo, z)} - E^^, (7) 

z 

with xq and yo lateral coordinates of the tip apex atom, and Ep^^ the Fermi energy of the 
tip material. 

In the above LDOS formulae, nT{E) and n'^{E) denote electron charge DOS projected 
to the tip apex and the ath surface atom, respectively, 

nriE) = nl{E) +n^p{E), 

n%{E) = nf{E)+nf{El (8) 

t and J, relative to their local spin quantization axes. Similarly, mp^E) and m'^{E) are 
electron magnetization DOS (MDOS) projected to the tip apex and the ath surface atom, 
respectively, 

mT{E) = nl{E) -n^(^), 

m-{E) = nf{E)-nf{E). (9) 

ifa is the angle between the spin moments of the tip apex and the ath surface atom. Above, 
the spin-resolved atom-projected DOS (PDOS) quantities, np^{E) and n'^^'^{E), are cal- 
culated from first principles. For this task any available ab initio electronic structure code 



can be used. Spin-resolved PDOS can also be calculated at finite temperatures, if we as- 
sume a Gaussian broadening of the peaks at the k-resolved spin-dependent electron energy 
(Kohn-Sham) eigenvalues, ^^^(fc), obtained at zero temperature, as 



= ^^_l_e-(^-4^-^te))VG^ f ^z^^MM^^^^M,i^^^ 

k 3 V 

^{tip apex atom} 

= ^^_l=e-(^-4'''(^))'/^' j rfVvl>f'^t(r)vI/f'^(r), (10) 



k j T/ 

^{ath. surface atom} 



with ^^-^■''(r) the spin-dependent electron wavefunctions corresponding to s^si!^) ^^P 
(T) and surface (S), respectively, and j the energy band index. The integral over the atomic 



volumes can be performed either in the atomic sphere or within the Bader volume [35|. In 
the present study we use integral over atomic spheres. The Gaussian parameter G could, in 
general, be temperature dependent. In our slab calculations, we fixed its value to 0.1 eV, 
which always provided smooth n^'^{E) and n'^'^'^{E) functions. 

LDOS can also be written in terms of energy dependent spin-polarizations, P{E), as 

LDOSixo, 2/0, z, E, V) = AEJ^ e-2"(^'^)l^"^("«'^°'")-^-lnr(E)ng(E) [l+PT{E)P^{E)cosip^] 

a 

(11) 

where the spin-polarization is defined as 

Pt{E) 



mriE) 


nl(E) - 


- n^^{E) 


TlTiE) 


nl{E) - 


f 4(E) ' 


m%{E) 


nf{E) 


-nf{E) 


n%E) 


nf{E) 


+ nf{E) 



PsiE) = = . ■ (12) 



Using Eq.(ll) of Ref. |29| and our LDOS expression, the differential conductance at the tip 



apex position and at energy E is 

^(xo, 2/0, z, E, V) = ji^Ef e-2«(^'^)li^"p(-0'^0'^)-^«lnT(E)n^(E) [1+Pt{E)P^{E)cos^, 

a 

(13) 

This is an important step to define dl/dV without the need for calculating the tunneling 
current. Multiplying the LDOS with AE results in a dimensionless quantity, which is mul- 
tiplied by the conductance quantum /h in order to arrive at our dl/dV expression. This 
means that nj'{E)AE electron states from tip and n'g{E)AE states from each surface atom 
contribute to the differential conductance at energy E, and in our model, dl/dV is propor- 
tional to the LDOS, which contains both surface and tip electronic information. Electronic 
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structure {PDOS) of tip apex can be calculated at the same level as surface electronic prop- 
erties, but it is also possible to combine different levels of electronic structure calculations, 
or include simplified model tip electronic structures into our approach. For example, assum- 
ing an electronically fiat maximally spin-polarized {Pt{E) = 1) ideal magnetic tip with e.g. 
nj-{E) = 1/AE, the differential conductance reads 

a 

Here, Eq.(jl]) has been assumed for the vacuum decay, and there is no V dependence. 

In single point STS experiments the tip is fixed above a surface atom. Its z position is 
determined by the bias set point V and the tunneling current /. Additionally, a modulating 
voltage Vm with a small amplitude is added to V, and dl/dV can be obtained from measuring 
current modulation Our model corresponds to this experimental setup, and the bias 
set point V determines the relative energetic position of sample and tip electronic structure, 
such that Ep^^ = Ep = Ep + eV , see Figure [1] for a sketch. The tip set point can be given 
either by V and J or by V and z because there is a one to one correspondence between z 
and / at given bias voltage, V , i.e. z{x, y) corresponds to the height profile of a constant 
current contour, / = Iq- Note that in all of our calculated single point dl/dV spectra we 
fix the tip apex z=3.5 A above the surface Fe atom. 

The general energy dependence of dl/dV in Eq.f lT3|l can be cast into bias dependence, 
U, which is varied, while V is fixed. Using E = Ep + eU we obtain 

^{xo,yo,z,U,V) = £!(AE)2^e-2'^(^^+'=^'^)l^"^(^«'^0'^)-^'^l (15) 

a 

X nriE^ + eU- eV)n%{E^p + eU)[l + Pt{EI + eU - eV)P^{E^ + eU)cosipa\. 

Note that by integrating our dl/dV expression with respect to eU in a given energy window 
defined by the bias voltage V and temperature, the tunneling current £aii^be calculated, 

26l |. except 



which expression is identical to the starting point of other STS theories js, 7, Is, 
the fact that our tunneling current is formulated within the atom-superposition approach. 
This will be reported in the future. In the present study, we focus on the simulation of 
SP-STS. Similarly to Eq.([T]), the differential conductance can also naturally be decomposed 
into a non-spin-polarized {TOPO) and a spin-polarized [MAGN) part, 

dl , - dlrOTAL , TT ^^TOPO , jj ... dluAGN , TT M\ 

^{xo,yo,z,U,V) = — — — {xo,yo,z,U,V) = {xo,yo, z,U,V)+ — — — {xo,yo, z,U,V), 



(16) 



thus, the contributions can be analyzed separately. This is also valid if 
we start from the tunneling current and define the dl/dV as its derivative 



with respect to bias voltage, see e.g. Ref. [26|, section III. A, where the oc- 
curring extra background term can be straightforwardly separated into TOPO 
and MAGN parts containing j dEn%E)-§^ [e-'^<^'^^\^Tip{^o,yo,z)-R^\^^^-^ _ ^y^^-^ 

cos(pa J dEmg{E)^ ^Q-ME,v)\RTipi^o,y(hz)-RJ^^(^]T; — eV^)] terms, respectively, inherent to 
our atom-superposition approach. 

The presented method for simulating STS can also be applied for nonmagnetic systems, 
where all magnetic contributions are equal to zero. Moreover, note that this method can 
be generalized in order to obtain two-dimensional differential conductance maps, however, 
without knowing the tunneling current, only in the so-called constant distance mode, where 
the tip apex position is varied in a plane parallel to the surface. This does not correspond 
to the usual experimental setup measuring tunneling spectra on a constant current contour, 
which ensures constant tip-sample distance even for rough surfaces [s], [ll| . The extension of 
our model in this direction is under way, and will be reported in a subsequent paper. 

III. RESULTS AND DISCUSSION 

We simulate differential conductance tunneling spectra of the Fe(OOl) surface, measured 
with an ideal, electronically flat, maximally spin-polarized tip, and with a model ferromag- 
netic Ni tip. The ideal tip has two advantages: (1) due to its electronically featureless 
character the bias set point does not play a role and the sole effect of tip magnetization 
direction on the spectra can be investigated, and (2) assuming maximal spin-polarization 
for the tip, the biggest magnetic effect on the spectra is expected. On the other hand, 
electronic structure of the Ni tip apex atom has been calculated from first principles. 

Tunneling spectrum of Fe(OOl) has been studied experimentally earlier using a nonmag- 
netic W tip 36|. According to that measurement the obtained peak at +0.17 V corresponds 
to a surface state of Fe(OOl) close to the T point. Ref. js^] reports a slightly different peak 
position of the same surface state at +0.3 V. This difference might be due to different tip 
geometry [i 



however, to the best of our knowledge, so far magnetic tips have not been 
considered for studying the Fe(OOl) surface. Here, we demonstrate that this surface state 
peak can be tuned either by changing tip magnetization direction or bias set point. The 
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results suggest that the sensitivity of STS depend crucially on the tip and it can be enhanced 
by finding favorable combination of these factors. 

Spin-polarized coUinear electronic structure calculations have been performed with stan- 
dard Density Functional Theory (DFT) methods within the Generalized Gradient Approxi- 
mation (GGA) implemented in the Vienna Ab-initio Simulation Package (VASP) 39|-|4l|. A 
plane wave basis set for electronic wavefunction expansion together with the projector aug- 
mented wave (PAW) method 42| has been applied, while the exchange- correlation functional 
is parametrized according to Perdew and Wang (PW91) 43|. 

We calculated a 9-layer Fe slab, where the surface layer on one side and the first subsurface 
layer have been fully relaxed. After relaxation the interlayer distances are reduced by 1.2% 
and 6.2% compared to bulk, respectively. We also checked a thicker slab of 13-layer Fe, where 
both surface layers have been relaxed and found no significant difference in the tunneling 
spectra. A separating vacuum region of 10 A width in the surface normal (z) direction has 
been set up between neighboring supercell slabs. For calculating the projected electron DOS 
onto the surface iron atom in our 1x1 surface unit cell we used anllxllxl Monkhorst- 



Pack (MP) [4^ k-point grid. The magnetic moment of the Fe surface atom is 2.84 /i^. The 
average electron workfunction of the iron surface is calculated to be (ps = 3.99 eV using 
Eq.(|6]). The spin quantization axis of Fe points along the direction of its in-plane magnetic 
moment. 

Figure [2] shows simulated differential tunneling spectra (solid lines) z=3.5 A above the 
surface Fe atom in the bias range from —1.0 V to -1-1.0 V using the ideal tip, employing 
Eq.(IT3j). The spectra are rescaled such that they can be shown together with the spin- 
polarization, Ps{E). The total tunneling spectra (TOTAL) are decomposed into topographic 
(TOPO) and magnetic [MAGN) contributions, according to Eq.f lT6|) . Focusing on the 
nonmagnetic contribution (TOPO, solid red curve), one can immediately observe that we 
obtain STS peaks at -1-0.2 V and slightly below -1-0.6 V, where the former one is in excellent 
agreement with previous experiment 36| , and both peaks are in good agreement with more 
sophisticated calculations employing a multiple scattering description for electron tunneling 



45|. These peaks originate from minority d electrons, in agreement with Refs. 



19 



Left 



and right parts of Figure [2] show spectra obtained with assumed tip magnetization direction 
antiparallel [C = cosip = —1) and parallel (C = cosip = +1) to that of Fe, respectively. Note 
that the total and topographic tunneling spectra are always positive in the whole energy 
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range, while the magnetic contribution can also be negative. Moreover, the topographic 
contribution equals in both parts as it is independent of tip magnetization, whereas the 
magnetic contribution changes sign, and the total dl/dV curve changes accordingly. While 
it is enhanced and its peaks are more pronounced in the antiparallel setup, it is lowered 
and considerably flattened in the parallel case, the peak at +0.2 V even disappears. This 
is an evidence that by modifying the magnetization direction of tip, the surface state peak 
can be tuned, and generally suggests that the sensitivity of STS on magnetic samples can 
be enhanced by choosing tip magnetization direction properly. For achieving this goal, and 
due to Figure [2], we identify the magnetic contribution, (IImagn / dV to play the key role. 

In order to understand this mechanism more, we study electronic properties of the surface 
Fe atom, Ps{E), {ns{E) + ms{E)cosip) , ns{E) , 'ms{E)cosip (dashed lines in Figure [2]). The 
calculated spin-polarization is —0.79 at the Fermi level, i.e. it is negative as has been reported 
recently by Ferriani et al. jl^, where the full-potential linearized augmented plane wave 
method has been applied. The shape of our energy dependent spin-polarization is also in 
good agreement with Ferriani's result. Comparing the remaining listed electronic properties 
to tunneling spectra it is clearly seen that they are highly correlated, i.e. dlTOTAi/dV oc 
{ns + niscosLp), dIxopo/dV oc ns, dluAGN / dV oc mscosip, each pair drawn with the same 
color in Figure [2l The multiplicative factor connecting each of these pairs is the energy 
dependent transmission coefficient for tunneling electrons, e"^''^^)^, where we applied Eq.(jl]) 
for the energy dependent vacuum decay. In Figure [2] the corresponding electronic properties 
are shown multiplied by e"^^'^-^^')^, and rescaled in the same way as dl/dV, such that e.g. 



and similarly for the rest of the correlated pairs. Note that the conversion between bias 
voltage and energy is ell = E — Ep, i.e. U = corresponds to Ep. Since the considered vac- 
uum decay, Eq.(jl]), is monotonously decreasing with increasing energy, the relation between 
the rescaled electronic properties and the corresponding dl/dV spectra is the following: For 
positive bias voltages. 



dluAGN 

dV 



(f/ = 0) = e 



ms {Ef) cosLf = e''-^^^^ms {E^) cosip, 



(17) 



dlroTAL 



(f/>0) > e 



2k (Sf 



> [ns [E^F + ef/) + ms {E^ + eU) COSif] , 



dV 




diMAGN 

dV 



(f/ > 0) > e 




ms [Ep + eU) coscp 



(18) 
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since e 2K(£;|+e!7)2 ^ ^ 2k(£;|)2_ Hg^g^ -y^^g employed Eqs. ffT^ and (fTSjl . Similarly, for negative 
bias voltages, e~^''(^^+''^)^ < e"^''^^^)'^, and, thus, 

dhoTAL^jj < 0) < e-2<^^)^ [ns + e?7) + ms + e?7) cos^] , 



dV 
dlropo 



(f/ < 0) < e"2"(^^)%5 (^1 + ef/) 



< e-2''(^^)"|m5(E| + e[/)cos(/?|, (19) 



This bias dependent relation between the correlated pairs can clearly be seen in Figure [2] 
(compare solid and dashed lines of the same color), which is, in turn, due to the energy 
dependent vacuum decay. 

From Figure [2] we identified the magnetic contribution, diMAGN / dV to be responsible for 
tuning the STS peaks. At a given energy E this is, in effect, the tunneling transmission 
coefficient (e"^'''^^^^) times the MDOS of surface Fe {ns{E)Ps{E)) times the cosine of the 
angle between spin quantization axes of surface and tip {cos(p). Since 6"^"^^)^ and ns{E) 
are always positive, the sign of diMAGN / dV is determined by Ps{E)cosip. Let us focus on 
the energy region E| - OAeV < E < E^ + l.OeV, where the spin-polarization is negative. 
If tip magnetization is antiparallel (parallel) to that of Fe, cosip = —1 [cosip = +1), then 
Ps{E)cosip is positive (negative), and so is dluAGN / dV . Adding this term to diTOPo/dV, 
which is not affected by tip magnetization, the total differential tunneling spectrum is en- 
hanced if tip magnetization is antiparallel, and decreased if it is parallel to Fe. This explains 
the main message of Figure [2l Although the surface state peak is already well obtained by 
using a nonmagnetic tip (dlTOPo/dV), the dl/dV signal can be further improved by setting 
the tip magnetization direction antiparallel to that of Fe. Note that in our ideal tip the 
energy independent spin-polarization was set to Pt = 1, thus, the largest effect on dl/dV 
occurs in this case. Similar, but reduced effect is expected if we set the tip spin-polarization 
to < Prp < 1, whereas dl/dV would be enhanced in the studied energy regime by applying 
a tip magnetization direction parallel to Fe if Pt < 0. Setting Pt = corresponds to a 
nonmagnetic tip, which results in dl / dV=dlTOPo/ dV . 

In Figure |2] we show simulated tunneling spectra considering contribution from one surface 
Fe atom only. However, in Eq.( fT4l) the summation over a should, in principle, be carried out 
over all surface atoms. Since the tunneling probability decays exponentially with increasing 
tip-sample distance, it is expected that a finite number of surface atoms should be enough 
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to be included in the summation in order to obtain converged dl/dV functions. In Figure 
[3] we study this convergence of the simulated differential tunneling spectrum above an Fe 
atom in the bias range from —1.0 V to +1.0 V using the ideal tip. We consider the case of 
antiparallel tip magnetization direction only, and the spectra are rescaled in the same way 
as in Figure [2J We show dl/dV spectra by including different number of surface Fe atoms 
in the summation over a in Eg. (1141) . i.e. one Fe atom in a 1 x 1 surface unit cell (IFe, dotted 
line), nine Fe atoms in a 3 x 3 surface cell (9Fe, dashed line), and twenty- five Fe atoms in 
a 5 X 5 surface cell (25Fe, solid line), where the spectrum is calculated above the central 
Fe atom with lateral coordinates (xo,yo) in each case. Since the magnetic surface unit cell 
is identical to the chemical unit cell, all surface Fe atoms have the same local electronic 
structure (PDOS). The spectrum obtained by one Fe contribution is the same as the one in 
the left part of Figure |2] drawn by black solid line. It is clearly seen that by including more 
atoms in the summation, the spectrum is growing with no change of the peak positions. We 
find that convergence is rapid, i.e. calculating dl/dV from a 5 x 5 surface cell is sufficiently 
converged, and contribution from all Fe atoms in a 7 x 7 cell means a relative increment 
of less than 10~^ compared to the 5x5 cell in the studied bias range. One should keep 
in mind that the independent orbital approximation for vacuum decay of electron states is 
employed, and taking into account orbital variations 3J| would alter the fine structure of 
our calculated spectra without changing the peak positions. 

We believe that within the model of the atom-superposition approach our finding to 
obtain convergence of the dl/ dV spectrum with respect to the spatial extension of the sample 
surface contributions is generally valid for any sample surface. The reason is the exponential 
factor describing electron tunneling transmission, which is decaying rapidly as the tip apex- 
surface atom distance increases. Considering heterogeneous sample surfaces, we always take 
into account the full chemical unit cell closest to the tip apex position for summation over a, 
and, thus, all important peaks appear in the dl/dV spectrum. Convergence of dl/dV with 
respect to spatial extension is obtained by involving atoms of neighboring full chemical unit 
cells in the summation over a. Moreover, note that in our Fe(OOl) surface by including more 
atoms in the summation does not change our findings for the effect of tip magnetization 
direction on the tunneling spectra described in Figure [2l 

In the above discussion we omitted energy variations of tip electronic structure for the 
purpose of studying only the effect of tip magnetization direction on the differential tunneling 
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spectrum. However, considering realistic tips the situation is somewhat more comphcated 
as diMAGN / dV at a given energy E is proportional to nT{E)PT{E)ns{E)Ps{E)cos(p. Let 
us study tunneling spectrum of the same Fe(OOl) surface by probing it with a ferromagnetic 
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25|. 



Ni tip. Such tips are routinely used in SP-STM and SP-STS experiments 

The Ni tip has been modeled by a 7-layer Ni film slab with (110) orientation, having one 
Ni apex atom on both surfaces, i.e. with a double vacuum boundary. According to previous 
findings it is sufficient to assume one tip apex atom on top of a metal surface as tip model, 
and there is no need to simulate more complex geometries, since the electronic structure of 



the tip apex does not change considerably, as it has been shown e.g. for an Fe tip [19|. For our 
purpose of demonstrating the effect of tip apex electronic structure on SP-STS spectra, such 
simple tip model is sufficient. In our tip the apex atom and the topmost surface layers have 
been relaxed on both sides. The interaction between apex atoms in neighboring supercells is 
minimized by choosing a 3 x 3 surface cell, and a 15.4 A wide separating vacuum region in z 
direction. Moreover, a 5 x 5 x 1 MP k-point grid has been chosen for obtaining the projected 
DOS onto the apex atom. Employing Eq.(I71), the local electron workfunction above the tip 
apex is 0t = 4.52 eV, and Eq.(l5]) has been used to determine the vacuum decay. 

Figure [U shows the calculated electronic structure, P{E),n{E),m{E) of the Ni tip apex 
atom (top part) and the Fe surface atom (bottom part) in the energy range [—2.5eV, +2.5eV] 
with respect to the corresponding Fermi energies. For Fe we find additional important 
PDOS peaks outside the [E^ - leV, E^ + leV] energy range reported in Figure O Peaks 
at Ep + 1.25eV and Ep + l.SeV have minority d character, while the peak at Ep — 1.5eV^ 
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m. 



originates mostly from majority d electrons, all these in good agreement with Refs. 
Focusing on Ni, we find that nx{E) is almost constant above Ep^^ + O.GeV, and most 
importantly, the spin-polarization is —0.91 at the Fermi level, Ep^^ , and it is negative and 
high in absolute value, i.e. |Pt(^)| > 0.8 between E]^^^ - 0.3eV and E^^^ + 0.3eV. The 
energetic relation of tip and sample electronic structures is determined by the bias voltage, 
V, illustrated in Figure[Il where ^=+1.0 V has been chosen. As Ep^^ = Ep + eV, the whole 
tip electronic structure is shifted by eV with respect to that of the sample. This means that 
depending on V, different electron states are involved in the tunneling process for calculating 
the tunneling current. Simulating differential tunneling spectra, the bias voltage is called 
bias set point and it fixes the relative energetic position of tip and sample electron states. As 
we learned from Figure [21 the total dl/dV signal can be tuned depending on the magnetic 
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contribution, (IImagn / dV , which is proportional to nT{E)PT{E)ns{E)Ps{E)cosip at a given 
energy E, in our model including energy variations of tip electronic structure. Depending on 
the bias set point the latter product can vary considerably, even can change sign, which, in 
effect, determines whether dl/dV is enhanced or decreased at the given energy. Since rixlE) 
and ns{E) are always positive, the decisive factor for the sign of the magnetic contribution 
is the effective spin-polarization, Pt{E)Ps{E)cos(p. 

In order to illustrate this effect, we consider three different bias set points taking con- 
tribution from one surface Fe atom to dl/dV in Figure HI This Figure shows simulated 
differential tunneling spectra (TOTAL, solid lines) and topographic contribution (TOPO, 
dashed lines) above the surface Fe atom in the bias range from —1.0 V to -1-1.0 V using 
our model Ni tip, employing Eq.f lT^ . depending on the bias set point, 0.0 V, -1-0.5 V, -1-1.0 
V, indicated by different colors. The spectra are rescaled such that they can be shown to- 
gether with the effective spin-polarization, Pnp{E)Ps{E)cosip, drawn by dotted lines. The 
topographic contributions for the corresponding bias set point are the same for both tip 
magnetization directions, and the iron surface state peak position is more or less reserved, 
while the peak heights change depending on the bias set point. This effect is due solely to 
bias set point, and not on tip magnetization, as dlxoPo/dV is proportional to nrp{E)ns{E) 
at energy E, and the bias set point determines the relative position of tip and sample elec- 
tronic structures. On the other hand, the total dl/dV curves change considerably depending 
not only on bias set point but also on tip magnetization direction, similarly to the obser- 
vation assuming an electronically fiat tip electronic structure, see Figure El Focusing on 
the antiparallel tip magnetization direction [C = cosip = —1, left part of Figure HI), it is 
seen that the iron surface state peak disappears at 0.0 V bias set point (orange solid line), 
it is shifted to —0.10 V by using -1-0.5 V (black solid line), while it remains at the same 
position and is enhanced by using -|-1.0 V (brown solid line). On the other hand, the effect 
is different by setting tip magnetization direction parallel to that of Fe (C = cos(p = +1, 
right part of Figure Hj). Here, for all considered bias set points the peak position remains, 
and for 0.0 V and -1-0.5 V the peak height is enhanced, while for -1-1.0 V the peak height is 
decreased. All these findings show an evidence that the STS peaks can be tuned by changing 
the bias set point. The reason for this effect is nT{E)ns{E) for the topographic part, and 
Pt{E)Ps{E)cos(p for the magnetic part, whose sign determines the sign of diMAGN / dV and, 
thus, the total dl/dV related to diTOPo/dV . This is clearly seen for all considered cases: if 



14 



Pt{E)Ps{E)cos(p > then dl/dV > dlroPo/dV in the corresponding energy regime, and 
similarly, if Pt{E)Ps{E)cos<^ < then dl/dV < dlroPo/dV. 

The results suggest that the sensitivity of STS on a magnetic sample can be enhanced by 
choosing the appropriate bias set point and magnetic tip. diTOPo/dV peaks are expected at 
local maxima of nT{E)ns{E), while the maximal diMAGN / dV can be obtained in case of hav- 
ing parallel tip and sample spin-polarization vectors with +1 value of the spin-polarization, 
each. This effect should also be observed in STS experiments. By designing the proper tip 



material, possibly wit 



study, or e.g. Ref. 



1 the help of electronic structure calculations, i.e. according to our 
19|, the advantage on the sensitivity of SP-STS measurements can be 



expected. However, the controlled preparation of magnetic tips for SP-STS experiments is 



not at all an easy task 



46|. Since our model assumes one tip apex atom, at present, we can 



not take rough tip structures or nanotips into account. On the other hand, it is naturally 
possible to simulate such rough tip structures within our atom-superposition framework by 
considering summation over different tip atoms contributing to dl/dV. This could be a 
research direction in the future. Though our present study does not help for improving tip 
preparation techniques, it is demonstrated that the governing factor of the fine structure 
and sensitivity of SP-STS spectra is the effective spin-polarization. Note that above con- 
siderations are generally valid for any combination of magnetic sample and tip, and could 
explain the observed difference in the bias dependent structure of the measured STS spectra 
above Co islands of opposite magnetization 24l. |25|. 

In the present study we omitted energy variations of tip and sample spin quan- 
tization axes. However, in certain combinations of tip and sample the situation 
is even more complicated as diMAGN / dV at a given energy E is proportional to 
nT{E)PT{E)ns{E)Ps{E)cos(p{E), i.e. the angle between spin quantization axes can also 
depend on energy. We will address this question in the future. 



IV. CONCLUSIONS 



Motivated by recent SP-STS experiments we studied tip effects on the tunneling spec- 
trum of a model magnetic Fe(OOl) surface. By considering an ideal electronically flat and 
maximally spin-polarized tip we found that STS peaks are sensitive to the tip magnetization 
direction relative to the surface. In case of a model Ni tip the role of the bias set point for 
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the tip is highlighted, which fixes the relative energetic position of sample and tip electronic 
structures. We showed evidence that the fine structure of the tunneling spectrum is governed 
by the effective spin-polarization. In conclusion, our results suggest that the sensitivity of 
STS on a magnetic sample can be tuned and even enhanced by choosing the combination of 
appropriate magnetic tip and bias set point. 
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FIG. 1: (Color online) Sketch of the relative energetic position of sample and tip electronic struc- 
tures. P(E),n{E),m{E) of the Ni tip apex and the surface Fe atom, and their relation at applied 
bias voltage, V are shown. Depending on V, different electron states are involved in the tun- 
neling process for calculating the tunneling current, while V (bias set point) fixes the relative 
energetic position of tip and sample electron states when simulating differential tunneling spectra, 
i.e. Ep^^ = AMPLE _j_ The energy scale is shown in the bottom left part (1 electronvolt), 
and according to that the bias voltage in the Figure is +1.0 V. 
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FIG. 2: (Color online) Simulated differential tunneling spectra (TOTAL) and their topographic 
(TOPO) and magnetic (MAGN) contributions (solid lines) z=3.5 A above an Fe atom on the 
Fe(OOl) surface, assuming a maximally spin-polarized and electronically flat tip using Eq. (|14|) 
and Eq. (|16p . Tip magnetization direction is antiparallel {C=cosip = — 1) (left part) and parallel 
{C=cosip = +1) (right part) to that of Fe. Note that in the two parts the topographic contribution 
is the same, while the magnetic contribution changes sign, and the total dl/dV curve changes 
accordingly. By modifying the magnetization direction of tip, the surface state peak can be tuned. 
For comparison, spin-polarization, Ps{E), and other electronic properties of the surface Fe atom, 
{ns{E) + ms{E) * cos(p),ns{E),ms{E) * cos(p, each multiplied by e~'^'^^^F)^ are shown in both 
parts (dashed lines). These properties correlate well with the corresponding spectra (same color), 
see text for details. 
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FIG. 3: Simulated differential tunneling spectra 3.5 A above an Fe atom on the Fe(OOl) surface, 
and its convergence depending on the contribution from the number of Fe atoms in the summation 
over a in Eq. (jl4p . i.e. from one Fe atom (IFe, dotted line) in a 1 x 1 surface unit cell (the spectrum 
is taken above this atom), from nine Fe atoms (9Fe, dashed line) in a 3 x 3 surface cell (the spectrum 
is above the central Fe), and from twenty- five Fe atoms (25Fe, solid line) in a 5 x 5 surface cell (the 
spectrum is above the central Fe). Note that each surface Fe atom has the same local electronic 
structure. For the simulation a maximally spin-polarized and electronically flat tip is applied using 
Eq. ()14p . and tip magnetization direction is antiparallel to that of Fe. 
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FIG. 4: (Color online) Simulated differential tunneling spectra (TOTAL, solid lines) and topo- 
graphic contribution (TOPO, dashed lines) 3.5 A above an Fe atom on the Fe(OOl) surface, using 
a model Ni tip (see text for details), and depending on the bias set point for the tip (curves with 
different colors). The chosen bias set points are given in parentheses. They fix the relative en- 
ergetic position of tip and sample electron states, see Figure [3l Tip magnetization direction is 
antiparallel {C=cosip = —1) (left part) and parallel {C=cos(p = +1) (right part) to that of Fe. 
Note that the topographic contributions for the corresponding bias set point are the same for both 
tip magnetization directions, while the total dl/dV curves change considerably, similarly to Figure 
[TJ This is due to PT{E)Ps{E)cosip variations depending on the bias set point, which are shown 
with dotted lines. Thus, by modifying the bias set point, the surface state peak can be tuned. 
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